UNCLASSIFIED 


ad.  291  803 


IZeptoduced- 
iuf  tUe 


ARMED  SERVICES  TECHNICAL  INFORMATION  AGENCY 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 

fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


6S~/~  6 


O' 


o 

00 

0> 

W 


MEMORANDUM 

RM-3148-PR 

DECEMBER  1962 


GO 

*g*'**’Vr* 

ca 

CD 

Lu 


s-J 


CD 

«C  v...;'" 

fC  C 

o  -• 


"■k* 


NOTE  ON 

A  HYPERSONIC  TAIL-SHOCK  PROBLEM 

R.  E.  Meyer 


PREPARED  FOR: 

UNITED  STATES  AIR  FORCE  PROJECT  RAND 


A  SIIS 


RAI1D 


\  .  •  ■ 


SANTA  MONICA  •  CALIFORNIA- 


;i  DEC-/ 

uL'vddl-,  u  w  uil^ 

TiS!A 


u 


MEMORANDUM 

RM-3148-PR 

DECEMBER  1962 


NOTE  ON 

A  HYPERSONIC  TAIL-SHOCK  PROBLEM 

R.  E.  Meyer 


uact  NoaAF  4qSr^Tr7nnhy  the  Unite1dlStatr  Air  Force  under  Project  RAND -Con- 
tract  Wo.  At  49(638 1-700  monitored  by  the  Directorate  of  Development  Planning 

Sd^n  thTS  Memora  ?,nd  Te^oSy,  Hq  USAF.  Views  or  conclusions  con 

nr  nnl-  *  f  fi.  i  r  ?  j  c  sh°llld  not  be  interpreted  as  representing  the  official  opinion 

of  this  Me  ^  L,nUed  State*  Alf  Force-  Permission  to  quote  from  or  reproduce  portions 
of  this  Memorandum  must  be  obtained  from  The  RAND  Corporation.  P 


-7 


R-pnD 


(?aifrox<xti4H 


Santa  monica  •  California  - 


t  700  MAIN  sr 


-iii- 


PREFACE 


Oliis  Memorandum  presents  a  preliminary  idealized  analysis  of  a 
part  of  the  flow  field  close  behind  a  hypersonic  projectile.  As  such, 
it  is  of  interest  primarily  to  those  engaged  in  hypersonic-wake  research. 
This  work  is  part  of  a  continuing  study  of  the  physical  ard  electro¬ 
magnetic  characteristics  of  hypersonic  trails  in  the  atmosphere. 

An  analysis,  by  the  same  author,  of  a  related  problem  was  pub¬ 
lished  as  RM-2929-PR,  Note  on  an  Axially  Symmetric  Expansion  Behind 


a  Body  in  Hypersonic  Flow. 
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SUMMARY 


For  the  purpose  of  studying  the  structure  of  the  inviscid  gas  flow 
near  the  narrow  neck  of  the  wake  behind  a  projectile,  the  formation  of 
a  shock  by  a  steady  gas  stream  converging  toward  an  axis  of  symmetry  is 
investigated.  The  flow  is  assumed  to  be  a  conical  field  with  vertex 
at  the  point  where  the  inner  boundary  of  the  stream  first  reaches  the 
axis  and  where  the  shock  starts. 

A  detailed  treatment  is  given  for  the  case  where  the  flow  down¬ 
stream  of  the  shock  is  uniform,  axial,  and  of  large  Mach  number  M,  so 
that  the  semi vertex  angle  r  of  the  shock  is  snail,  but  Mr  is  not.  It 
is  also  assumed  that  the  shock  is  at  least  fairly  strong.  An  additional 
small  parameter  then  appears,  and  a  uniform  first-order  approxin»tion 
with  respect  to  it  is  derived.  The  streamlines  are  found  to  be  nearly 
straight  and  parallel,  but  the  gas  is  found  to  suffer  a  precompression, 
upstream  of  the  shock,  such  that  the  inner  boundary  of  the  inviscid 
gas  stream  is  a  vacuum  line.  Similarly,  vacuum  is  approached  asymptot¬ 
ically  on  all  streamlines  with  distance  upstream  from  the  shock. 

It  is  shown  that  the  result  concerning  the  occurrence  of  vacuum 
must  be  expected  to  be  independent  of  all  but  one  assumption--that  the 
gas  impinges  on  the  axis  at  a  non-zero  angle.  Accordingly,  tail-shock 
formation  in  two-dimensional  and  axially  symmetrical  flows  must  be  ex¬ 
pected  to  differ  in  significant,  still-unknown  ways. 
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I.  INTRODUCTION* ** 


When  a  body  of  revolution  with  a  blunt  base  flies  at  constant, 
supersonic  speed,  the  wake  is  known  to  form  a  neck  from  which  a  tail- 
shock  issues  (Fig.  l). 


Fig.  1 

If  the  boundary  layer  is  thin,  the  neck  is  very  narrow,  especially  in 
the  axially  symmetrical  case,  and  the  shock  forms  at  a  radius  that  is 
small  compared  with  both  the  radius  of  the  base  and  the  radius  of  curva 
ture  of  the  shock.  The  flow  near  the  shock  and  just  outside  the  region 
in  which  viscosity  plays  a  strong,  direct  role  my  therefore  possess  a 
structure  similar  to  that  obtained  in  the  limiting  case  of  zero  neck 
and  boundary -layer  thickness.  In  this  case,  the  shock  would  be 


* 

The  author  is  indebted  to  Prof.  J.  D.  Cole  for  drawing  his  at¬ 
tention,  after  completion  of  the  work  reported  here,  to  a  note  by 
A.  A.  Nikol'skii  on  the  boat- tail  problem.  Since  this  note  is  un¬ 
translated  and  not  easily  accessible,  its  main  argument  is  reproduced 
in  Sec.  VI.  Nikol’skii's  argument  complements  the  work  reported  here 
so  as  to  leave  no  escape  from  the  conclusion  that  axially  symmetrical/ 
tail-shock  formation  is  not  well  understood. 

** 

At  high  Mach  numbers,  the  wake  does  not  look  thin,  but  if  Prandtl' s 
approximation  concerning  the  momentum  balance  in  the  wake  is  relevant 

then  the  limit  of  infinite  Reynolds  number  should  remain  similarly  rele¬ 
vant. 
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expected  to  form  on  the  axis  at  a  point  embedded  in  a  region  in  which 
the  radius  takes  values  that  are  large  compared  with  the  mean  free  path 
but  only  small  compared  with  any  other  lengths  representative  of  the 
flow.  Accordingly,  the  existence  of  a  continuum  flow  with  conical  sym¬ 
metry  (Fig.  2)  should  be  expected  in  this  region.  This  flow  will  be 
studied  in  the  following. 


Fig.  2 

It  might  be  anticipated  that  a  significant  difference  exists,  just 
as  in  the  analogous  nose-shock  problem,  between  the  cases  of  two-di¬ 
mensional  and  axially  symmetrical  flow.  In  the  former  case,  the  so¬ 
lution  would  be  expected  to  consist  simply  of  a  plane  shock  flanked  by 
two  regions  of  uniform  flow.  In  the  axially  symmetrical  case,  on  the 
other  hand,  a  nontrivival  conical  solution  may  exist  such  that  the 
shock  deflection  angle  t5  (Fig.  2)  differs  from  the  approach  angle  T0. 

A  distinction  then  arises  between  the  ratio  p^/p^  the  limiting  pres¬ 
sures  on  the  two  sides  of  the  shock  and  the  ratio  p  /p  of  the  pressure 

__  c  EL 

Pg  just  downstream  of  the  shock  to  the  pressure  p&  on  the  boundary  OA 
(Fig.  2)  of  the  incident  stream.  These  are  respectively  the  shock- 
pressure  and  the  approach -pres sure  ratios.  In  order  to  determine  the 
latter  in  terms  of  the  approach  angle,  an  explicit  solution  will  first 
be  derived  for  a  restricted  range  of  circumstances.  In  Sec.  V,  it 
will  then  be  shown  that  the  main  results  are  independent  of  most  of 
these  restrictions. 
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II.  GOVgRNING  EQUATIONS 


For  the  quantitative  treatment,  we  assume  that  the  flow  downstream 
of  the  shock  is  purely  axial.  Since  the  conical  assumption  indies 
uniformity  of  the  entropy  except  for  a  jump  across  the  shock,  it  follows 
that  the  flow  downstream  of  the  shock  is  uniform  For  hypersonic  flight 
speeds,  the  shock  angle  x  may  be  taken  to  be  a  small  parameter.  It 
will  also  be  assured  that  the  Mach  number  M  downstream  of  the  shock  is 
large,  so  that  Mr  is  not  small.  The  problem  then  falls  within  the 

hypersonic  small- disturbance  theory. 

Accordingly,  let  x,  r,  i,  p,  u,  and  v  denote  respectively  the  axial 
distance,  radius,  pressure,  density,  and  axial  and  radial  velocity  cam- 
pone  nts,  and 


x  =  x,  r  =  rr,  u  =  u^(l  +  x2u),  v  =  xu^vCr, x) 
P  =  P27>^T2p(r,  x),  p  =  P2p(r,  x) 


-where  the  subscript  2  refers  to  the  uniform  flow  downstream  of  the  shock. 
Then  as  x  ->  0,  the  first  approximations  to  the  equations  of  continuity 

v.  (2) 

and  momentum  become 


(2) 


5v  .  v9v  ,  1  op 
3x  +  3r“  +  p  ^ 


0 


(3) 


The  former  can  be  satisfied  by  the  introduction  of  a  stream  function 
i|f  such  that 

=  r-p,  ^  =  -rPv  (4) 


and,  since  viscous  effects  are  neglected 


“  cn(^)  (5) 

P7 

where  y  is  the  ratio  of  the  specific  heats  and  is  assumed  to  be  constant 
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Let  denote  the  Mach  number  just  upstream 
let  a  =  sin2  (t  +  tS),  where  t&  is  the  stream 
shock  (Fig.  2).  The  density  and  pressure  ratios 


of  the  shock,  and 
deflection  across  the 
across  the  shock  are 


P 

P 


2 

1 


and  for  a  perfect  gas 


(3) 


l_ 

Ra 


Pa 


R 


2  +  (y  -  l)a 

"T7  +  lja 


p  _  gycr  -  (7+1) 

7  +  1 


(6) 


let  the  subscript  1  denote  conditions  just  upstream  of  the  shock.  Then, 
by  Eq.  (l) 


Px  =  R, 


“L  = 


1 

7M2t2PR7 


(7) 


The  continuity  of  normal  mass-flow  rate  and  tangential  velocity  across 
the  shock  implies  (Fig.  2) 


tan  r  =  R  tan  (t  +  t6) 


and  so  to  the  first  order 


v 


1 


=  -8 


(8) 


The  final  boundary  condition  required  is  the  continuity  of  ^  across 
the  shock. 

The  assumption  of  conical  symmetry  implies  that  ^  has  the  form* 
t  =  |  r2F(n),  n  =  f  (9) 


•nr  P 

To  avoid  confusion,  note  that  i|/  =»  x  f(r/x)  in  Ref.  1  (Eq.  (13b)) 
is  misprinted.  Since  the  line  x  =  0  falls  in  the  nonuniform  flow 
region  of  the  tail- shock  problem,  the  use  of  7)  is  more  convenient 
here. 
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and  by  Eqs.  (3)  and  (5),  F  must  satisfy 


(10) 


upstream  of  the  shock,  where  ui  =  const  =  ,  and  the  boundary  condi¬ 

tions  at  the  shock  are 


F( 1 )  =  1,  F’(l)  =  2(1  -  R) 


(11) 


by  Eqs.  (l),  (8),  and  (9).  The  mathematical  problem  is  to  integrate 
Eq.  (10)  from  T|  =  1  towards  smaller  values  of  T|  in  order  to  find  the 
value  T1  =  -0  for  which  F  =  0,  and  to  compute  the  pressure.  The  anal¬ 
ogous  problem  for  the  nose  shock  was  solved  numerically  in  Ref.  1.  For 
the  tail  shock,  however,  an  analytical  approach  will  be  found  more 
suitable.  Since  it  is  the  interval  of  F,  not  of  T),  which  is  known  a 
priori 


will  be  regarded  as  the  independent  variable,  and 

G(F)  =  R  +  (l  -  R)T1  (12) 


as  the  dependent  variable.  The  system  of  Eqs.  (10)  and  (ll)  then 
transforms  into 


F2G"  +  I  FG'  = 


/^(G  -  R) 

(1  -  R)2 


G  -  K' 


[<°  - R)0" +  c'2]  [F  -  -ion 


7+1 


(13) 


G(l)  =  1,  G'(l)  =  | 


(lb) 


where  primes  denote  differentiation  with  respect  to  F. 
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III.  STRONG- SHOCK  APPROXIMATION 

The  solution  of  Egs.  (13)  and  (l4)  is  seen  to  depend  on  the  two 
parameters  R(cr)  and  0^(0).  At  hypersonic  flight  speeds,  the  experi¬ 
mental  evidence  indicates  tail  shocks  which  possess  appreciable  pressure 
ratios  near  the  axis.  The  shock  strength  0  -  1  =  sin2(x  +  x8)  -  1 
might  therefore  be  expected  to  be  large  for  a  significant  range  of  flight 
conditions.  The  limit  of  R  as  l/o  -*  0  depends  on  the  chemical  state 
of  the  gas  and  might  be  small,  but  for  a  perfect  gas  with  7  =  1.4, 

R(»)  *  1/6;  accordingly,  R(»)  will  not  be  regarded  here  as  a  «m~n 
parameter.  On  the  other  hand,  P  1  ->  0  as  l/a  — »  0,  and  so  0,  corre¬ 

sponding  to  the  well-known  fact  that  the  entropy  upstream  of  a  shock 
becomes  negligible  compared  with  that  downstream,  in  the  limit  of  in¬ 
finite  shock  strength.  Moreover,  the  relevant  parameter  will  not  be 
but 

(l  -  R)2  ^t2P(1  -  R)2 

and,  for  example,  if  7  =1.4  (1.2)  and  shock  strength  0  -  1  =»  2  only, 
corresponding  to  a  shock  pressure  ratio  of  3.33  (3.2),  the  value  of 
7^t  p.  is  only  0.27  (0.19).  The  assumption  that  u  «  1  should  there¬ 
fore  be  proper  for  an  important  range  of  hypersonic  flight  conditions. 

Hence,  let  ^  -> 0.  Then  by  Eg.  (3),  v  =  const  on  the  streamlines 
as  well  as  on  the  conical  rays,  so  v  =  const  upstream  of  the  shock, 
just  as  in  the  two-dimensional  case,  and  6  =  8  *»  (l  -  R)/R.  However, 
this  approximation  yields  no  information  on  the  approach  pressure  be¬ 
yond  5a  «  i2.  By  Egs.  (4),  (8),  and  (9),  moreover,  p  =  RE1/2,  indi¬ 
cating  vacuum  at  the  approach  line  AO  (Fig.  2). 

Now  consider  the  problem  of  obtaining  a  higher  approximation.  The 
first  approximation,  corresponding  to  tUj  =  0  and  v  =  const,  is 

G^T)  ~  E1/2  (l6) 

If  the  perturbation  problem  were  a  regular  one,  the  second  approxi nation. 
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G2(F),  would  be  obtained  from  Eqs.  (13)  and  (14)  by  substituting  G1 
on  the  right-hand  side  of  Eq.  (13);  i.e.,  solving 


+  I  °2  "  t  "  R)?(7//2)"2 


(17) 


with  G0( 1)  =  1,  G^(l)  -  1/2,  which  yields 

)0 

(i .  A*)-1)] 


G  (F)  =  F1/2  +  — H — 
2  2-7 


[(Hf R-1 


R 


F1/2  )  +  ^  1  -  F(7-1 

(18) 


3-7 


But  as  F  -»  0,  the  right-hand  side  of  Eq.  (17)  tends  to  ®  if  n  ji  0,  and 
hence,  while  Eq.  (18)  provides  an  accurate  description  of  the  flow  in 
most  of  the  region  between  AO  and  OS  (Fig.  2),  there  yet  remains  a 
narrow  "boundary"  sector  adjacent  to  the  approach  line  AO  in  which 
Eq.  (16)  fails  to  furnish  even  a  first  approximation  for  arbitrarily 
small  n  >  0. 

The  root  of  the  difficulty  is  easily  understandable  from  Eqs.  (3) 
and  (5),  which,  if  q  and  s  denote  respectively  the  velocity  magnitude 
and  arc  length  on  the  streamlines,  may  be  combined  to 

3v 

q  31  pdy 

The  nondime  ns i o nal  pressure  gradient  represented  by  the  right-hand 
side  is  singular  at  0  (Fig.  3)  but  small  everywhere  else,  so  v  changes 
very  little  if  any  chosen  streamline  is  followed  upstream  from  the  shock 
over  a  finite  distance  (Fig.  3). 


r 


Fig.  3 
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2 

But  to  reach  very  small  values  of  the  conical  variable  F  =  2\[r/r  re¬ 
quires  following  the  streamline  over  a  very  large  distance,  over  which 
the  effect  of  the  small  pressure  gradient  might  accumulate  significantly. 
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IV.  THE  BOUNDARY  SECTOR 

The  results  of  the  preceding  section  show  that  as  long  as 
F  =  2ijr/r2  is  not  small,  the  first  approximation  is  given  by  Eq.  (13) 
with  the  right-hand  side  deleted,  so  that  G  =  G1(F)  =  F1/2.  But  in 
the  boundary  sector  where  F  and  hence  also  G  are  sufficiently  small, 

Eq.  (13)  must  reduce  to  a  different,  approximate  equation  in  which  u 
plays  an  essential  role. 

Note  that  this  sector  near  the  boundary  F  =  0  must  contain  a  tran¬ 
sition  sector  for  which  Eq.  (16),  while  inaccurate,  yet  provides  the 
correct  relative  orders  of  magnitude.  In  this  transition  sector,  there¬ 
fore,  |F|  «  1,  | G |  «  R,  and 

| FG ' j  «  R  (19) 

so  that 

P  “  F  “  ^2^  ~  |  ,  H  =  2G'  (20) 

and  Eq.  (13)  must  reduce  to 


(F2H7+1  -  hR)H'  -  -  |  H^FH7  +  u)  (21) 

The  transition  sector  begins  where  F  is  small  enough  for  u  to  be  signi¬ 
ficant  in  Eq.  (21);  by  Eq.  (19),  that  must  be  when  F2H7+1  =  0(u),  while 

•V 

still  FH'  »  u.  Hence,  the  limiting  form  of  Eq.  (13)  must  become 

(F2H7+1  -  uR)H’  -  -  |  FH7+2  (22) 

whence 


4o 

HT  = 


4yRR 


5-7 


(5  -  7) 


+  const 


(23) 


and  const  =  1  in  order  that  Eq.  (23)  be  identical  with  Eq.  (l 6)  in 
the  "regular"  sector.  Hence 


F2  = 


4uR 


3  -  7H 


7+1 


(2k) 


-lo¬ 


in  the  transition  sector.  The  regular  sector,  then,  corresponds  to 


F  >>  h  « 


and  the  transition  sector  corresponds  to 

7  =  o(h-2)  -  o(^A3-r)) 


(25) 


In  the  range  expressed  in  Eq.  (25)>  FH  =  ±.e.f 

Eq.  (19)  remains  valid  for  1  <  7  <  3.  By  Eq.  (22),  moreover 

H<r+1)/2  i  - ■  o<«1/2>  ^-y),z 


0(^(7-! )/( 6-2 7)) 


(26) 


in  the  range  of  Eq.  (25),  and  so  the  corresponding  range  of  G  is  «  R, 
and  hence  still  |g|  «  R.  Note  also  that  fV+1  decreases  as  7  does, 
by  Eqs.  (24)  and  (2 6),  while  H  increases.  But  Eq.  (24)  shows  that 
fV+1  possesses  a  lower  bound  in  the  transition  sector 

jV+1  2  Ju£L 

3-7 

and  so  in  Eq.  (21),  F^H^"*”1  -  p.R  £  (7  +  l)uR/(3  -  7).  Hence,  however 

p  /(  ^-7' )  y 

small  the  fraction  of  p.  '  ' J  '  to  which  F  may  have  decreased,  FH  »  p 

and  Eq.  (22)  remains  the  first  approximation  to  Eq.  (13) J  Accordingly, 

Eq.  (24)  is  the  uniform  first-order  approximation  to  the  solution 

dG/dF  =  H/2  of  Eqs.  (13)  and  (14)  for  1  a  F  a  0. 

The  sector  of  Eq.  (25)  thus  describes  the  transition  from  the 

"regular"  sector,  where  Eq.  (l 6)  is  the  first  approximation,  to  the 

2 /(3-y) 

boundary  sector,  F  «  p  '  ' J  ,  where  the  first  approximation  is 

fV+1  -  Jtt£- 

3-7 


(27) 


It  follows  from  Eq.  (2 6)  that  G(0)  «  1,  and  the  first  approximation 
to  the  approach  angle  remains 


The  first  approximations  for  the  nondimens ional  density  and  pressure 
in  the  boundary  sector  F  «  ^/(3~7)  are 


On  the  approach  line  AO  (Fig.  2),  the  density  and  pressure  vanish,  not 
because  the  shock  strength  is  infinite,  but  because  the  conical  solu¬ 
tion  implies  a  precompression  along  the  streamlines  from  vacuum  at  in¬ 
finity  upstream,  to  the  density  p1  and  pressure  p^  at  the  shock. 
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V.  A  QUALITATIVE  ARGUMENT 


It  may  be  instructive  to  review  the  problem  solved  in  the  pre¬ 
ceding  sections  in  the  light  of  a  simpler,  but  only  qualitative,  argu¬ 
ment  which  illuminates  especially  the  role  played  by  the  assumption 
of  conical  symmetry.  This  assumption  states  that  p  and  v  are  to  be 
functions  of  only  a  single  independent,  conical  variable  and  implies 
that  the  entropy  measure  oi  is  constant  upstream  of  the  shock.  In  the 
following,  the  independent  variable  will  be  chosen  to  be,  not  J\  =  x/r, 
but  the  arc  length  a  on  an  arbitrarily  chosen  streamline  measured  up¬ 
stream  from  the  shock  in  multiples  of  some  suitable  reference  length. 
Then,  if  q  =  (l  +  v^)  represents  the  velocity  magnitude 

$£.  |£  Tic  d£ 

ox  1  -  vT]  do  *  dr  1  -  vT)  do 

Hie  same  relations  hold  for  the  derivatives  of  v.  The  momentum  and 
energy  equations,  Eqs.  (3)  and  (5),  may  therefore  be  combined  to  obtain 


dv 


7-2  _ IL 


1  -  vTl 


do 


(28) 


and  the  continuity  equation,  Eq.  (2),  may  be  written  with  the  help  of 
Eq.  (28)  as 


dp  v 
p  do  ~  rD 


(29) 


I)  =  1  -  ^-1  - £ _ 

(1  -  vTl)2 


(30) 


Equation  (29)  may  be  used  in  turn  to  cast  Eq.  (28)  into  the  form 


3. 

v 


dv 

dp 


=  7“P 


7-1 


— h _  i_ 

1  -  vT]  rD 


(31) 


Observe  that  1  -  vTl  =  0  on  the  approach  line  AO  (Fig.  2)  and  hence 
that  l-vT]-»Oasa— >coon  our  streamline,  far  upstream  of  the  shock. 
Similarly,  l/r  ->  0  as  a  ». 
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We  now  introduce  the  strong-shock  assumption  o>'«  1,  with  R  not 
small.  Since  T)  =  1  and  v  =  -(l  -  R  )/r  at  the  shock,  by  Eq.  (8),  Eq. 
(31)  shows  v(a)  to  be  approximately  constant  as  long  as  1  -  vt)  =  0(1). 
To  this  approximation,  therefore,  the  streamline  is  straight,  and 
q/(l  -  vt))  is  a  linearly  decreasing  function  of  c r,  while  r  is  a  line¬ 
arly  increasing  function  of  a;  therefore,  Eq.  (29)  shows  p  to  decrease 
in  inverse  proportion  to  some  positive  power  of  r  as  the  streamline  is 
followed  upstream  from  the  shock. 

Sufficiently  far  upstream,  however,  we  must  anticipate  that  a  will 
become  so  large  that 

1ct2  =  0(1) 


This  relation  defines  a  "transition"  segment  of  the  streamline,  on 
which  cop  a  «  1  still,  so  that  v  remains  approximately  constant,  by 
Eq.  (31  )>  while  the  law  of  dependence  of  p  on  0  changes  appreciably, 
by  Eq.  (29).  When  this  new  trend  only  begins  to  take  effect,  D  must 
still  be  positive,  and  since  v  <  0,  p  must  still  decrease  with  increas¬ 
ing  a. 

As  the  streamline  is  followed  farther  upstream,  there  are  there¬ 
fore  just  two  possibilities.  One  is  that  a  stage  is  reached  where 
y- 1  2 

cop  0  »  1;  that  must  certainly  occur  if  p  is  bounded  away  from  zero 

on  the  streamline.  But  before  this  stage  can  be  reached,  a  value  of 
0  must  occur  at  which  D  =  0  and  dp/du  and  dv/d 0  are  singular,  by  Eqs. 
(29)  and  (31).  Since  0  represents  also  a  conical  variable,  the  singu¬ 
larity  is  a  limit  line,  and  the  solution  is  multivalued.  A  solution 

of  this  kind,  ^ut  with  expansion  in  the  stream  direction,  is  envisaged 
by  Reyn . '  ' 

The  second  alternative  is  that  D  remains  positive,  and  hence 
p  -*  0  as  0  -*  co  upstream. 

It  should  be  observed  that  the  argument  just  given  is  independent 
of  the  assumption  of  uniform,  axial  flow  downstream  of  the  shock.  What¬ 
ever  that  flow,  if  the  assumption  of  conical  symmetry  is  made,  the  shock 
must  be  straight,  and  if  its  inclination  T  «  1  (Fig.  2)  and  1^T  is  not 


The  quantitative  discussion  of  Sec.  Ill  confirms  that  pr_1/(l  -  vn)2 
increases  with  0  when  1  -  vr)  =  0(l).  u 


small ,  where  is  the  Mach  number  just  downstream  of  the  shock,  then 

Eqs.  (2),  (3),  and  (5)  are  valid,  and  hence  also  Eqs.  (2 9)  and  (3l)* 

The  argument  then  remains  quite  unchanged,  unless  the  highly  implausible 
assumption  be  made  that  the  gas  enters  the  shock  at  a  glancing  angle 
«  t  (so  that  1  -  v1  «  l).  Again,  since  only  the  product  dp'  appears 
on  the  right-hand  side  of  Eqs.  (29)  and  (31),  the  argument  is  unchanged 
should  the  density  ratio  R  =  pl/^2  ”  P1  across  tlie  be  small,  even 

though  v  /  0(1)  in  this  case.  In  fact,  the  argument  is  based  essentially 
on  the  structure  of  the  denominator  D;  therefore,  it  is  also  independent 
of  the  strong-shock  assumption.  The  value  of  D  at  the  shock  is 


by  Eqs.  (5)  and  (l),  and  if  the  plausible  assumption  be  made  that  the 
gas  is  still  converging  towards  the  axis  when  it  enters  the  shock  (so 
that  v1  <  0),  then  the  condition  ^  1  is  sufficient  to  insure  >  0. 
It  follows  from  Eq.  (29)  that  a  compression  precedes  the  shock,  and  if 
v  remains  negative  as  the  streamlines  are  followed  upstream  a  limit  line 
D  >=  0  can  be  avoided  only  ifp  ->0as  (l  -  vl])  -»  0  far  upstream. 
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VI.  NUCOL  *  SKIT 8  S  ARGUMENT 

A  quite  different  argument^  shows  that  the  hypersonic  assumption 
is  also  unnecessary  and  that  not  even  the  conical  assumption  is  essen¬ 
tial.  The  characteristic  equations  of  steady,  axially  symmetrical  flow 
of  inviscid  gas  are 


+  cot  adp  —  sin  0  sin  u  dr 
q2p  sin  (e  +  4)  r 


on  -  tan  (0  +  4)  (32) 


where  0  is  the  stream  direction,  4  is  the  Mach  angle,  and  where  ordi¬ 
nary,  not  hypersonically  distorted,  variables  are  used.  More  explicitly, 
for  a  perfect  gas  with  constant  specific  heat,  Eq.  (32)  may  be  written^5 * 


d0  ±  dn>  ±  ~  dS  +  sin  0  sin  4  cosec  (0+4)  —  =  0 


on 


dr 

dx 


tan  (0+4) 


(33) 


where  7(7  -  1)\  =»  sin  4  cos  4,  S  is  the  specific  entropy,  and 

»/2 

is  the  Prandtl  angle,  which  is  a  monotonic,  continuously  differentiable 
function  of  4  with  <n(0)  -  (n  -  1)jc/2,  n2  -  (7  +  l)/(7  -  1). 

Consider  now  the  flow  past  a  body  of  revolution  with  boat  tail 
(Pig.  4),  the  family  of  Mach  lines  C'E'  sloping  down  to  the  body  surface, 
and,  in  particular,  the  limiting  member  CE  of  this  family  which  meets 
the  body  surface  at  the  tip  of  the  tail.  If  the  flow  is  not  totally 
degenerate,  the  respective  differences  between  the  values  of  0,  0,  and 
S  at  the  points  C  and  E  must  be  finite.  On  the  other  hand  if  Eq.  (32) 
or  (33)  is  integrated  along  the  Mach  line  from  C  to  E,  the  term  pro¬ 
portional  to  r  furnishes  a  divergent  contribution,  unless  0  -  0  or 
P  ■  0  at  E.  But  the  meridian  slope  of  the  body  is  negative  at  E  if  the 
"tail  is  not  cusped. 
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